Survival of quantum effects for observables after decoherence 
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When a quantum nonlinear system is linearly coupled to an infinite bath of harmonic oscillators, 
quantum coherence of the system is lost on a decoherence time-scale td- Nevertheless, quantum 
effects for observables may still survive environment-induced decoherence, and be observed for times 
much larger than the decoherence time-scale. In particular, we show that the Ehrenfest time, which 
characterizes a departure of quantum dynamics for observables from the corresponding classical 
dynamics, can be observed for a quasi-classical nonlinear oscillator for times t ^ td- We discuss 
this observation in relation to recent experiments on quantum nonlinear systems in the quasi-classical 
region of parameters. 

PACS numbers: 



In the last few decades there has been extensive theo- 
retical, and more recently experimental, research on the 
quantum-classical transition. It has been noted that ev- 
ery physical system is, in fact, an open quantum sys- 
tem interacting with its environment. Consequently, the 
evolution of the reduced density matrix of the system 
(obtained after tracing over the environmental variables) 
evolves in such a way that quantum coherent effects are 
quickly suppressed. This process of environment-induced 
decoherence has been considered to be an essential in- 
gredient of the quantum-classical transition On the 
other hand, despite the huge number of papers on this 
subject, only few of them deal with quantum nonlinear 
systems (e.g. iSSiQ). 

We consider in this paper the dynamics of a quantum 
nonlinear oscillator (QNO) 



(1) 



interacting with a bath of linear oscillators which are ini- 
tially in thermal equilibrium. Here a (a^) are annihilation 
(creation) bosonic operators, u is the linear frequency 
and fjL is the parameter of nonlinearity. The QNO is ini- 
tially prepared in a coherent state in the quasi-classical 
region of parameters. In the classical limit (a — > a. 



0, \a\ 



J - an action of the 



classical linear oscillator) the Hamiltonian becomes 
Hci = wJ + /iJ^. In what follows we use the following 
dimensionless notation: r = ut, fl — hfi/uj, fici = fJ-J/cu 
and e = h/J. Thus, the quantum parameter p. can be 
presented as a product of two parameters, quantum and 
classical, p, = e/id. The parameter /id characterizes the 
nonlinearity in the classical system, and can be written as 
/^ci — {J/'2(^){du!ci/dJ), where Wci = dHd/dJ — lu + 2/iJ 
is the classical frequency of nonlinear oscillations. The 
limit /ici <C 1 corresponds to weak nonlinearity, while 
/ici > 1 corresponds to strong nonlinearity. The parame- 
ter e is a quasi-classical parameter. Namely, e ~ 1 corre- 
sponds to the "pure" quantum system, and e <C 1 corre- 



sponds to the quasi-classical limit, which is the subject 
of this paper. 

We study the following problem: What are the pa- 
rameter conditions for observation of quantum effects on 
expectation values (observables) in the QNO dynamics. 
We describe the dynamics for the QNO for observables 
taking into account five characteristic time-scales which 
naturally appear in this system. Three of them char- 
acterize the time-scales of the QNO evolving under the 
Hamiltonian dynamics: (i) Td = 27r/a;d - the period 
of nonlinear classical oscillations; (ii) te - the so-called 
Ehrenfest time, which indicates the characteristic time- 
scale at which quantum dynamics for observables starts 
to depart from the corresponding classical dynamics; (iii) 
Tji - a quantum recurrence time, which describes the 
time-scale for quantum recurrences of observables under 
the Hamiltonian evolution. There are also two charac- 
teristic time-scales related to the interaction of the QNO 
with the thermal bath: (iv) r^i - a decoherence time, 
which characterizes the decay of the non-diagonal ma- 
trix elements of the reduced density matrix in the eigen- 
basis of the non-interacting Hamiltonian; and (v) Tj - 
a time-scale of relaxation of quantum observables due 
to the interaction with the thermal bath. We demon- 
strate that even if the decoherence time is much smaller 
than the Ehrenfest time, td te, one can still observe 
quantum effects for observables. Actually, the impor- 
tant condition for observation of quantum effects related 
to the Ehrenfest time-scale is te < Tj, which may be 
realized in modern experiments in the quasi-classical re- 
gion of parameters (e <C 1). This means that gener- 
ally the environment-induced decoherence is insufficient 
for recovering the quantum-classical correspondence for 
observables in quantum nonlinear systems. This is an 
important observation for at least two reasons: (a) It 
means that pure quantum effects can be observed for 
times much longer than to and (b) Pure quantum dy- 
namical effects can be important in experiments even in 
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the quasi-classical region of parameters. Finally, the clas- 
sical limit appears in our system under very natural con- 
ditions: Td ^ Tcl <^ <^ te Tr. 

One type of systems that could be considered for 
observing quantum effects in the quasi-classical region 
of parameters are Bose- Einstein condensates (BECs). 
These are particularly suited to analyzing the interplay 
between nonlinear dynamics and environmental interac- 
tions in the realm of quantum mesoscopic systems. This 
is because they are macroscopic matter waves, often de- 
scribed theoretically by the Gross-Pitaevskii (GP) equa- 
tion, which formally is a classical nonlinear field-theory. 
Going beyond GP allows one to understand the role of 
quantum effects in the quasi-classical region of param- 
eters. Another important feature of BECs is that they 
are experimentally easily accessible and controllable by 
means of trapping potentials and tunable interactions. 
BECs have already been used to demonstrate nonlinear 
dynamics of collapses and revivals of a coherent matter 
wave in the pure quantum regime Here, the conden- 
sate was initially trapped in the lowest energy band of a 
three-dimensional optical lattice. By adiabatically rais- 
ing the heights of the barriers it was possible to suppress 
tunneling between sites and at the same time maintain 
the system in the superfluid regime, so that within each 
lattice site independent coherent states were engineered 
with an average number of atoms of order one. Every 
BEC in each lattice site is well described in the single- 
mode approximation by the QNO Hamiltonian given in 
Eq. 1^, with oj being the trapping frequency in each lat- 
tice site and ji being proportional to the s-wave two-body 
scattering length. Other systems that could be used for 
observation of quantum nonlinear effects in the quasi- 
classical region of parameters are micro- 9] and nanome- 
chanical Iff] resonators, high-frequency cantilevers [Tlj . 
nonlinear optical systems and superconductive junctions. 

The QNO is one of the simplest quantum nonlinear sys- 
tems for which the breakdown of the quantum-classical 
correspondence can be exactly calculated. The quantum 
and classical dynamics of an initial coherent wave packet 
evolving under Eq. |Q were computed by Berman et al 
and by Milburn 3] . The characteristic time-scale for 
departure of the quantum dynamics from the correspond- 
ing classical one, the so-called Ehrenfest time, was intro- 
duced for this system in Q . In Q a similar problem was 
studied using the Q quasi-probability distribution. It was 
shown that the presence of non-positive definite second- 
order terms in the quantum evolution equation for Q, not 
present in the evolution of the classical probability den- 
sity, are responsible for quantum recurrences and prevent 
the appearance of fine-scale-structure "whorls" predicted 
in the classical description. In the interaction of the 
QNO with the environment (modeled by a thermal bath 
of harmonic oscillators) was studied in the limit of small 
nonlinearity, and it was argued that such interaction was 
effective in destroying quantum interference effects and 
restoring the classical phase-space structure. However, 
as we already stated, environment-induced decoherence 



is in fact ineffective in recovering the quantum-classical 
correspondence for this nonlinear system. 

In the following we compute quantum observables 
in the coherent state basis. For an arbitrary opera- 
tor function / = /(a^,a), the time-dependent expecta- 
tion value of such a function (observable), f{a*,a,t) — 



-iHt/n 



a), for an initial coherent state, Jo) 



satisfies a partial differential equation of the form [12| : 
df /dr — Kf, where K = Kd + tiKq. Here the operator 
Kci includes only the first order derivatives and describes 
the corresponding classical limit, while the other opera- 
tor Kq includes higher order derivatives and is responsi- 
ble for quantum effects. For the model given by Eq. 
we get 



9/ 
dr 



2fl\a\^) 
9(a*)2 



a 



d 
da* 

do? 



d_ 

da 



f + 



/• 



(2) 



In particular, for f — a the evolution of /(r) corre- 
sponds to the evolution of the condensate matter-wave 
field a{T) = (aja(r)|a). In this case Eq. ^ can be 
solved exactly : 



a(r) 



exp[|Q;p(e' 



!)]■ 



(3) 



The quantum evolution of this expectation value de- 
parts from the corresponding classical evolution adir) = 
ae-"^-'^ as a^r) = aci(r)e-^'/2r| _,_ _^ 
Odap/I'^r'^)], where te is the Ehrenfest time-scale given 
by 
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TE = 



2fl\a\ 



(4) 



The amplitudes of quantum and classical observables co- 
incide at times multiple of the quantum recurrence time- 
scale 



(5) 



Note that the quasi-classical limit, which is considered 
in this paper, corresponds to the following condition: 
te/tr, — £^/^/27r ^ 1. So, in what follows we will 
be interested only in the region of parameters where 
Te Tji. Quantum recurrences of the matter wave field 
of a BEC in the pure quantum regime a ~ 0{V) (or 
£ w f ) in each lattice site were observed in 7] at ~ 100 
{tR — 0.55ms), larger than the corresponding Ehrenfest 
time-scale te ~ 15. At the same time, quantum dynami- 
cal effects in the quasi-classical region of parameters have 
still not been observed in BECs. 

Expressing a = -^/j/^e^'^, we can rewrite Eq. as 



^ = (l + 2/ieO^ 



2e/ic 



d^f 
dJdB' 



(6) 



The quantum term appears as a singular perturbation 
of the classical equation because the small parameter 



£ multiplies the higher order derivative. Note that the 
quantum effects for observables vanish in two cases: (i) 
e = 0, which corresponds to the classical limit, and (ii) 
/ici — 0, which corresponds to the quantum linear oscil- 
lator. The fact that for nonlinear quantum systems the 
terms with high order derivatives in the evolution equa- 
tions for the density matrix and for the Wigner function 
represent a singular perturbation to the classical limit 
(Liouville function) is well-known. However, in spite of a 
large number of papers on this subject, from this fact it 
is still unclear what are the conditions for the quantum- 
classical correspondence for observables. The solution (jS)) 
of Eq. ^ for the observable a{T) (and also for an arbi- 
trary observable demonstrates that quantum effects 
(second order derivatives in Eqs. (|2l6|l ) represent a sin- 
gular perturbation to the classical equation for observ- 
ables, which includes only the first order derivatives and 
can be solved by the method of classical characteristics 
1^, 12] . This results in a secular behavior of quantum cor- 
rections in the solution for the observable ^(t) So, 
the question is: Under what conditions does the envi- 
ronment "kill" (if at all) the quantum corrections which 
represent a singular perturbation to the observables of 
the classical world? 

Following 0,13 we model the environment as a bath of 
harmonic oscillators in thermal equilibrium at a rescaled 
temperature (3 — TiujIk^T linearly coupled through posi- 
tion to the QNO [13. In the Born-Markov approxima- 
tion, the master equation for the reduced density matrix 
reads 



dp 



Ft,,,{p)+F^{p)+F,{p), 



where the first term, 

Firccip) = -i[d^a + fl{a^a)'^,p], 



(7) 



(8) 



corresponds to the free, unitary evolution, the second one 

Fvip) = ^[a + a^{il(r)a + aUl(T) + 

^{A2{T)a^d^A2{T)),p}], (9) 
accounts for dissipation, and the third one, 

i{MT)a-a^B2{T)),p]], (10) 

is related to noise. The time-dependent, operator- valued 
coefficients Ai and Bi depend on the frequency operator 

Cl^ l + pi(l-f 2ata), 

and on the spectral density of the environment 
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FIG. 1: Comparison between the numerically obtained deco- 
herence time td and the approximate relation given by Eq, 
I19II . obtained retaining only one term in the master equation 
The dashed line corresponds to r_D (theory) = rij (numerics) 
Data are obtained by varying parameters in the following 
regions: 
10"^ < 7 < 10 



20 < Jo < 100, 10"^ < M < 4, 10"^ < P < I 



that we chose to be Ohmic, with A a UV cut-off and 7 a 
system-environment coupling constant. Explicitly 



Mr) - 


/ dsri{s) cos{Cls)] 
Jo 


(11) 


Mr) = 


/ dsr]{s) sm{0,s); 
Jo 


(12) 


Mr) = 


/ dsi^{s) cos{fls); 
Jo 


(13) 


Mr) - 


/ dsvis) sm{ris), 
Jo 


(14) 



where the dissipation and noise kernels are respectively 
given by 



77(3) = 



Hs) 



dLO — J[ijj) sin(a;s); 

TT 



(15) 



du> — J [oj) coth{ — ) cos{cjs). (16) 

TT 2 



The matrix elements of the operators Ai and Bi can be 
straightforwardly computed in the Fock basis and shown 
to have an analogous behavior to that of the fl = 0, 
quantum Brownian motion case |l4j . 

To study the decoherence effects of the environ- 
ment we start by considering an initial Schrodinger cat 
state formed by large-amphtude coherent states, p(0) 



AA(|a) + |/?))((a| + (/3|), with 



> 1, and J\f a 



normalization constant. For simplicity we take a and /3 
to lie along a common radius, and we parameterize them 
as a = xe^^ and (3 = {x + 6x)e^^ . Decoherence is due to 
the term in the master equation containing Bi. In the 



coherent state basis, the off-diagonal matrix elements of 
the density matrix evolve as 



dT 



{aim « -2B,{Io,T){Sxf{a\m, (17) 



where Bi{Io,t) {a\Bi{T)\a) with Iq = \a^\. For r > 
1/A, Bi{Iq,t) is approximately equal to its asymptotic 
value 

i?i(/o,oo) = -_^coth(^-j, 

with = 1 + jx{l + 2/o). Therefore, the decoherence 
time-scale is 



1 



TD 



2B'^{oo){5xf 



(18) 



This decoherence time-scale coincides with the time-scale 
for exponential decay of quantum recurrences for an ini- 
tial coherent state \a) coupled to the thermal bath. In- 



deed, in this case {5x) 



In, so that 



TD 



tanh(/317/2) 



(19) 



We checked this estimation by numerical simulations of 
Eq. O using as initial state a quasi-classical coherent 
state. As can be seen in Fig. ^ the agreement is fairly 
good within the numerical errors. In the limit of small 
temperature /3fl ^ 1 and small non-linearity p,Io ^ 1 
the decoherence time-scale coincides with the one de- 
rived in 0; however, we stress that in the general case, 
differently from 0, the decoherence time, as given by 
(|19|l . depends on the parameter of nonlinearity /i. In the 
high-temperature limit, 13^1 <^ 1, the decoherence time 
is: td ~ Tjh^Lu/AkBTJci, where — 2/7 is the time 
scale of the relaxation of quantum observables due to the 
interaction with the environment. 

Having the density matrix elements, we can easily de- 
termine the average values of any observable, for in- 
stance, for the position x{t) = {x{t)) = (a(T) -|- 
a*(r))/-\/2- An example of such simulations is given 
in Fig. 12 The time-scale for the overall decay of the 
amplitude of recurrences, shown in Fig. is set by 
the decoherence time-scale T]j: the relative heights of 
two peaks, taken at two neighbor recurrent times, is re- 
duced by a factor exp(— ta/tc). An enlargement of the 
first bump of Fig. |3l is given in Fig. ^p. As one can 
see in both cases reported in the figure, the time-scale 
which governs the envelope of x(r) is the Ehrenfest time 
te, which is independent of the coupling to the environ- 
ment. The same is true for any of the following revival 
bumps, see Fig. Let us also notice that in Fig. 
the curves for 7 — 10^^ (solid line) and for 7 = 10^^ 
(dashed-dotted line) are slightly shifted one to the other. 
This is due to the fact that the frequency of motion 
is slightly dependent on the bath-oscillator interaction 
strength, = fl^ - ^K^/(h? + Vi^). This is not at all 




FIG. 2: The average position as a function of the dimen- 
sionless time r. In all cases Jo = 50, /3 = 1. a) Parameters 
are fi = 0.1 and 7 = 10~*. The dashed curve corresponds to 
exp(— t/td), where td = 18 > te ~ 0.7. b) is an enlarge- 
ment of the first bump of a). Two more curves have been 
added: a dashed-dotted line corresponds to the average po- 
sition for 7 — 10~^, so that td = 0.18 < te; and a dotted 
line that corresponds to the envelope exp(— r^/2r|;). c) is an 
enlargement of the third bump of a). The dotted line cor- 
responds to the envelope exp(— Tij/ro) exp[— (r — tr)'^ /2te], 
where tr — lOvr. d) Parameters are fi — j = 10~^, so that 

TD <. TE <. < TR. 



surprising since the renormalization of the frequency is a 
feature of the considered master equation [l^l • 

In Fig. 2d the average position is plotted for a case in 
which r£i ^ <C < {td — 0.74, Td = 3.12, te — 
7.1, T-y = 200 and tr = 314). This figure represents the 
most important result of the present work: Usually, in the 
quasi-classical region of parameters, and for rather large 
values of the coupling to the environment (7 > 10~^), 
the characteristic decoherence time-scale is much shorter 
than the Ehrenfest time scale, te, <^ te- Despite this, 
the system does not become entirely "classical", since 
quantum effects persist up to the Ehrenfest time. Indeed, 
for a;(0) = (a -f a*)/V2 = ^/2h = 10 and td = 0.74, the 
dependence a;(0) exp(— t/td) would give us, for example, 
for T = 10 the value 1.3 x 10~^, which is significantly 
smaller than the correspondent value 3.7 defined by the 
function a;(0) exp(-rV2r|) for x(0) = 10, r 10 and 
= 7.1 (which corresponds to the results presented in 
Fig. 2d.) 

In our model the classical limit corresponds to the fol- 
lowing inequalities: td ^ Td <^ t^ te tr (see Fig. 
|3Jl. Because in the quasi-classical region of parameters 
the inequalities te ^ tr and td <C Tc\ are always satis- 
fied, the really important condition for the classical limit 
is Tj <^ Te- In this case the system effectively behaves as 
a classical damped oscillator, and quantum effects cannot 
be observed. For comparison we plot in Fig. |31the overall 
decay of oscillation, given by r-y , with the one that would 
be given by the Ehrenfest time. The perfect agreement of 
the decay relaxation time with the data and their wrong 
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20 40 60 80 100 120 

T 

FIG. 3: The average position x(r) in the "classical" limit: 
td <€. Tci <^ ^ te <^ TR. Parameters are: p, = 10"'*, 
/3 = 1, 7 = 0.01, Jo = 50, so that td = 0.92, r^i = 2n, 
T-y = 200, TE = 707, TR ^ n xW^. The correct decay (dotted 
line), as given by the relaxation time r-y, is compared with the 
"wrong" decay (dashed line) given by the Ehrenfest time te- 



dependence on the Ehrenfest time is a manifestation of 
the classicality for this case. 

To gain a qualitative understanding of these numerical 
findings we consider a simplified version of our master 
equation Eq. Q at zero temperature, in which we keep 
only the effect of dissipation and decoherence due to the 
environment, and make the rotating wave approximation. 
In this way we get the standard master equation in quan- 
tum optics for the QNO: 

~ — i[a^a+/i(ata)^, p] + -^(2a/5at — a^ap— pa^a). (20) 

UT 2 

This equation is precisely the one considered in |j| , where 
an exact solution for the quasi-probability Q function 
was obtained, assuming an initial coherent state \a). In 
particular, an exact expression for time-evolution of the 
average position {x{t)) ~ ((a(T)) + (a'''(T)))/V2 can be 
written, where 

(a(T)) = ae'-^+^^'^e e (21) 

with k — 7/2/2. Note that for no coupling to the envi- 
ronment (7 = 0) we recover Eq. From Eq. H21() 
we can read the decay factor of the average position 
x{t) oc e~-^^'^\ It is given by 

^W-f + ^fftx (22) 
2 4^^ + 7^ 

(1 - e-'''' cos(2/Ir)) - ^e"'^^ sin(2/lT) . 
2^ 

Let us first analyze the case 7/2 ^ p,, corresponding 
to Te <C Ty. Let us express the time r around a given 
recurrence time as r = titji+t, where n is a non-negative 
integer. Assuming that the time r is much shorter than 




(O/CO , co/co 



FIG. 4: The Fourier spectrum of the average position versus 
the rescaled frequency, where ujd = 1-1- 2/i/o is the classical 
frequency. As the dashed black line we show the theoretical 
expression Eq. 12511 for the closed system. Parameters /3, fi, Iq 
are the same as in Fig. ©, while a) 7 = 10"^ td = 180 > 
TE = 0.7, b) 7 = 10-*, TD = 18 > TE ^ 0.7 c) 7 = 10~^ , 

TD = 1.8 > TE = 0.7 d) 7 = 10-^ TD = 0.18 <TE = 0.7 As 

one can see, decreasing of the decoherence time has no effects 
on the width of the Fourier spectrum. 



the relaxation time, <^ 1, and that /2f ^ 1, we can 
expand D{t) « |Q;p(2/i^f^ -|- nTR'y). Hence, in these 
limits, the decay of the average position is 

x{t) CX e-^'/2r| ^-nm/ro^ (23) 

where te is defined in Eq. Q and td — l/^\a\'^ is the 
zero temperature limit of Eq. H19I) for <C 1. There- 

fore, the decay of x{t) within any recurrence bump is 
determined by the Ehrenfest time-scale, and even in the 
limit Te ^ te the decay within the first bump (n = 0) is 
still governed by the Ehrenfest time-scale, which agrees 
with our numerical results presented above. This implies 
that some quantum effects survive the loss of quantum 
coherence due to the interaction with the environment. 
One the other hand, when 7/2 ^ Ji (i.e., t^ <C te) the 
classical limit is attained, and the decay is governed by 
the relaxation rate. 

Persistence of quantum effects after the decoherence 
time can be also observed analyzing the Fourier spectrum 
of the average position in time x{t) = x^e^'^'^. When 
the system is closed, i.e. no coupling to the environment, 
the Fourier components are given by 

a / uj - 1\ 

where 

A(n) = ^ r^%(e^-^-l)l"l^-2.Anr^^^ (24) 

7^ Jo 
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Estimating this quantity in the hmit /iT ^ 1 we obtain 
the frequency spectrum for the QNO, 



1 



: exp 



{ll> - Ucl)'^ 



x(0), 



(25) 



which is a Gaussian distribution centered around the clas- 
sical oscillation frequency = 1 + 2/2|ap with a spec- 
tral width given by the inverse of the Ehrenfest time- 
scale, Auj — Tj^^. In Fig. 2] the Fourier spectrum for 
the QNO coupled to the environment is shown for the 
case te <^ t^. As one can see, in all cases presented in 
Fig. 21 (whatever the relation is between the decoherence 
time and the Ehrenfest time) the width of the Fourier 
spectrum is always given by the inverse Ehrenfest time. 
When Tj <C Te the width of the spectrum is given by the 
relaxation rate, and the classical limit is obtained. 

The important condition for survival of the quantum 
effects for observables related to the Ehrenfest time-scale 
is <C Ty, which can be written in the form 



6 ^ ^ = 2McieV2 

TE 



T^ > 1. 



(26) 



For BECs the parameter of nonlinearity /id = iiJ/uj can 
be written as /ici = N y^aFmIoj2TTh, where N is the 
number of particles in the condensate, a is the s-wave 
scattering length, m is the mass of the atoms and to is 
the trapping frequency. The quasi-classical parameter is 
£ — 1/N. Therefore 
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BEC 



2mujN 



irfi 



-T^ > 1. 



(27) 



For example, for a = 5nm, m = 1.5 x 10~^^kg, uj/2tt = 
lOOHz and estimating the dimensionless relaxation time 
T-y from the lifetime of the condensate (say = Isec, so 
=ujt^ =27rxl02), we need a total number of particles 
^ 1. In the case of a cantilever (or a mechanical res- 
onator) the quasi-classical parameter is e = where n 
in the average number of levels involved in the coherent 
state of the cantilever. For the dimensionless relaxation 
time we take t^ = 2Q, where Q is the cantilever qual- 
ity factor. Then, for a cantilever the condition (|2(j|l takes 
the form 



e, 



cantilever 



> 1. 



(28) 



We take the following dimensional parameters the 
amplitude of the cantilever oscillations = lOnm, the 
spring constant kc — 6 x 10~^N/m and the frequency of 
the fundamental mode of the cantilever uJc/2tt — 6.6kHz. 
In this case, the number of cantilever levels can be es- 
timated as n w kcxI^/TiLOc ~ 6 x 10^^. We also take 
Q = 10^. Then, we have from Eq. the estimate 

for the parameter of nolinearity: /id ^ 0.2. We hope 
that these conditions can be experimentally realized, and 
quantum effects related to the Ehrenfest time-scale can 
be observed in the quasi-classical region of parameters. 

We acknowledge useful discussions with G. Nardelli. 
This work was supported by the Department of Energy 
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Defense Advanced Research Projects Agency (DARPA). 
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